A statistical portrait of the entanglement decay of two-qubit memories 
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We present a novel approach to the study of entanglement decay, which focuses on collective 
properties. As an example, we investigate the entanglement decay of a two-qubit system, produced 
by local identical reservoirs acting on the qubits, for three experimentally and theoretically relevant 
cases. We study the probability distributions of disentanglement times, a quantity independent 
of the measure used to quantify entanglement, and the time-dependent probability distribution of 
concurrence. Analytical results are obtained for initially uniformly distributed pure states. The 
calculation of these probability distributions gives a complete insight on how different decoherence 
channels affect the entanglement initially contained in the set of two-qubit pure states. Numerical 
results are reported for randomly distributed initial mixed states. Although the paper focuses 
in Markovian noisy channels, we show that our results also describe non-autonomous and non- 
Mar kovian channels. 

PACS numbers: 03.65.Ud,03.65.Yz,03.67.Mn 



I. INTRODUCTION. 

Entanglement is a key resource for quantum informa- 
tion processing and communication, which can be ma- 
nipulated, broadcasted, controlled and distributed 
A better understanding of the effects of decoherence on 
entanglement could accelerate the development of tech- 
nological applications of quantum information protocols. 
Its study under realistic, dissipativc conditions, is of fun- 
damental importance to assess the resilience of quantum 
information processing. The dynamics of quantum coher- 
ence and entanglement can be strikingly different. Un- 
der the influence of a noisy environment, coherences usu- 
ally decay asymptotically in time, while entanglement 
can do so in a finite time 0, Q . Extensive work [j 
including several experimental demonstrations [12 



has been devoted to this remarkable phenomenon, some- 
times called entanglement sudden death (ESD). A com- 
plex, yet partial picture of the entanglement dynamics, 
has emerged. Recently, a complete characterization of 
the dynamics of a two-qubit system in which only one 
of the qubits is coupled to a noisy channel [l^ has been 
presented and general approaches, based on geometrical 
arguments [l6j and quantum trajectories 17 1, have been 
proposed. 

In this work, adopting a broader perspective, we render 
a comprehensive description of the entanglement dynam- 
ics in two-qubit systems when corrupted by environment- 
induced decoherence. Our approach relies on a rather 
novel tool in the area: the statistical distributions of con- 
currence and of the ESD times. Obtaining analytic ex- 
pressions for them is far from being a symple task. In 
fact, previous analytic work was restricted to averages 
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jTsf and often relied on approximations [T^. Other re- 
ported explorations of entanglement evolution in bipar- 
tite and multipartite systems [20j . similar in spirit 
to our approach, are both numerical. Our analytic for- 
mulation, employed before in a kinematical context [2l|, 
allows us to depict a highly detailed image of the entan- 
glement time evolution. 

This paper is organized as follows. The dynamical evo- 
lution of a two-qubit system coupled to local channels 
is introduced in section |lll where the probability den- 
sities of the disentanglement times and of concurrence 
are defined. We consider three models of decoherence 
for a single qubit and assume identical environments for 
each qubit. In the next four sections the results for the 
probability densities of the disentanglement times and of 
concurrence are evaluated for the chosen environments, 
assuming initial pure states uniformly distributed. The 
results for mixed states, a single noisy channel, and non- 
autonomous systems are briefly presented in sections lVIIi 
IVIIII and IIXI respectively. Section |X] is devoted to the 
demonstration that, although we concentrate on Marko- 
vian processes with constant decay rates, our results also 
apply to non-Markovian processes. Some conclusions are 
drawn in the last section of this paper. 



II. DISENTANGLEMENT PROBABILITY 
DISTRIBUTIONS 



We investigate the global entanglement dynamics of 
two-qubit systems coupled to local reservoirs (chan- 
nels), a situation likely to hold for spatially separated 
qubits. The qubits are assumed to be quantum mem- 
ories; they do not follow independent unitary dynam- 
ics nor do they interact with each other. When- 
ever the initial system-environment state is separable 
with vanishing quantum discord (2^ [23j . the dynam- 
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ical evolution of the state of the two-qubit system is 
given by a completely positive map (24 . [25| . Then, 
p{t), the system state at the physical time t, is given 

by = E^,{4'^®Ef)p{Q)(^E^'^^<E>Ef^) = 

At(p(0)), where the time-dependent Kraus operators 
^(«) „ E\°'\t), acting on the a-th qubit, satisfy the 

trace-preserving condition X^i ^i"^^^!'^^ = a = 1,2. 
These operators can be experimentally determined by a 
quantum process tomography (26l | . 

We chose pairs of identical environments to act on the 
qubits, and regard the experimentally relevant cases of 
depolarizing (D), amplitude-damping (AD), and phase- 
damping (PD) channels. These generally non-Markovian 
channels are valid for finite or infinite environments and 
weak or strong coupling [27l |. The Kraus operators for 
these noisy environments, in terms of the (reparameter- 
ized) time g, are (0 < q < 1): En[q) = ^Jl - 3g/4I 
and Ei{q) = ^qj^Oi, i = 1,2,3, for case D; Eq ~ 
|0)(0| + Vr ~^|l )(l| and El = Vg|0)(l| for case AD; 
and So = VT^I, El = (0|, and E^ = (1| 

for case PD. Here cri,(T2, and 173 denote, as usual, the 
three Pauli matrices. Unless otherwise stated, we focus 
on the Markovian scenario, with constant decay rates, 
when q{t) = 1 — exp(— 7f), 7 being a positive constant. 
The motivation behind our selection of channels lies in 
their qualitative different long-time behavior, g — > 1, de- 
termining the dynamical evolution of entanglement on a 
gross level [iBl • All states will suffer ESD if the asymp- 
totic state belongs to the interior of the set of separable 
states S (case D). Some states will separate asymptoti- 
cally if the asymptotic state belongs to dS, the border 
between S and E, the set of entangled states (case AD). 
The case PD displays a set of stationary separable states, 
diagonal in the standard basis, which includes 4 pure 
states in dS. 

Rather than compute the concurrence (or any other 
entanglement measure) for specific initial conditions, we 
pursue a statistical approach to entanglement dynamics 
in open systems. Although randomly chosen normalized 
pure states = ^/-oo |00) + V'oi |01) + V'lo |10) -f |11> 
are used for our calculations, other choices might be 
sensible. We choose pure states because they gener- 
ally offer higher quantum correlations than mixed states, 
a key point for quantum information protocols. More- 
over, they are a convenient idealization of the almost 
pure states that can be produced under experimental 
conditions. If minimal prior knowledge about the state 
of the system is assumed [2^, or invariance under uni- 
tary transformations is required, it becomes natural to 
use the uniform (Haar) measure for all possible pure 
states of the system. Hence, the probability to find a 
state in a small volume (Pip = Yiij dipij dTp*j around ip is 
p{iP)d'^i) = T{A)5{l-{iP\^) )/7r4d2^, where T{x) denotes, 
as usual, the Gamma function. 

We characterize the disentanglement process through 



the reparameterized separation time probability density. 



Piqs)= d^ijd{qs-qsiip))pW, (1) 



where qs{ip) is the ESD time, i.e., the time at which 
p{q) = AgdV') (V'l), becomes separable. This distribu- 
tion is independent of the measure used to quantify the 
entanglement of the system. However, p{qs) may overes- 
timate entanglement if many states sustain small values 
of entanglement for long times, before undergoing ESD. 
Therefore, a complete characterization of entanglement 
decay also requires its quantification. 

The global evolution of entanglement is characterized 
by using the concurrence probability density. 



p{C; q) = / 5{C - C{q; V'))p(^) , (2) 



where C{q; ip) is the Woottcr's concurrence [2^ 
of the state p{q) = Aq{\ip) (-01). Concurrence is 
given by C{q;ip) = max{0, ^/A7 - ^^1' 
where are the eigenvalues of the matrix 

p{q) (4^^ ® ^) p{qy (4^^ ® 4^^) , being Ai the 

largest one. Here p* is the complex conjugate of p in the 
standard basis, and (T2 is the second of Pauli matrices. 
Then, for our chosen ensemble of initial states, p{C; q)dC 
is the probability to find a state with concurrence value 
in an interval dC around C at time q. In particular, at 
the initial time 9 = 0, eq. ([2|) provides the concurrence 
distribution over an ensemble of uniformly distributed 
pure states, p(C; 0) = 3CV1 - [HI, [30|. The con- 
tribution of the states with vanishing concurrence to 
the distribution p{C;q) can be singled out by writing 
p{C;q) = p(C;g) + S{q)S{C), where p(C; q) is the 
probability density for strictly positive concurrence and 
S{q) = 1 — /p " dC p{C;q) is the probability to find a 
separable state at a time q, where Cm = CM{q) is the 
maximum concurrence at that time. 

We will show that p{C; q) and p{qs) provide a com- 
plete portrait of the entanglement behavior for all times, 
enabling us to monitor entanglement degradation, con- 
centration, and even loss. Despite the dissimilar dynam- 
ics of the D, AD, and PD channels, in all cases we find 
analytical expressions for the disentanglement time and 
the concurrence, for initial pure states, and use them to 
evaluate the probability distributions of separation times 
p{qs) and of concurrence p(C; q) as a function of the repa- 
rameterized time. Our results were checked using purely 
numerical methods. 
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III. PROBABILITY DISTRIBUTIONS OF THE 
DISENTANGLEMENT TIMES. 

The disentanglement times, 

= 1 - l/v/l + 2Co, (3) 
4^°)=Co/2|^nP, (4) 

9r^-l-(v/Co' + rf^-Co) /d (5) 

for identical depolarizing, amplitude-damping and phase- 
damping channels, respectively, are found from the con- 
dition of separability, C{qs,ip) = 0. Here, we have 
set d = 4|'0oo'0oiV'ioV'iir^^, and Co = 2|?/>ooV'ii - 
'0oi'0io| = C(0, and have used superscripts indicat- 
ing the case/channel we are considering. We use Eq. ([T]) 
and these expressions for the separation time to obtain 
the following probability distributions of disentanglement 
time, 

(D). . _ 3gs(2 - qs)y/{q's - ^qs + 2) (3g^ - 6qs -f2) 




Dimensionless Separation Time qs 

FIG. 1. (Color online) Disentanglement-time probability dis- 
tributions; solid line (black) for local depolarizing channels 
(D); dashed line (blue), including a delta function contribu- 
tion, at gs = 1, for local amplitude damping channels (AD); 
and dotted line (red) for local phase damping channels (PD). 



IV. DEPOLARIZING CHANNELS. 



^(AD) 



arctan((75) 2-1-- 



2qs{l + qir 



2g| 



-Siqs-l), 



P 



(PD) 



(95) 



d6 ds dr Pj{s, r)S{qs 



In the latter case the dependency on the an- 
gle 9 comes from the initial concurrence Cq = 
V^y^s^ + r2 + (s2 _ 7.2) cos(e'). Here, Pj {s,r) ^ 

^^==^K ^ Iztr'-^ ^ stands for the joint probability dis- 
tribution of r and .s, with s + r = 2|'0oo'0ii| and s — r = 
2|V'oiV'io|- The complete elliptic integral of the first kind 
and elliptic modulus k is denoted by K{k^). 

The probability distributions of the separation times 
(Fig. [T]) display completely different behaviors. All 
present a rather conspicuous peak at around qs — 0.38, 
qs = 1, and qs = 0.59 for cases D, AD, and PD, respec- 
tively. However, the first is the widest, and the second 
is the more narrow and is actually a delta distribution. 
While in the first case all states become separable before 
qs = 1- 1/a/3 « 0.42, only ^ « 35.73% of the states 
experience ESD in the second case, and a zero-measure 
set of states separate asymptotically in the third case. 
The three cases we have considered can be tell apart by 
their disentanglement time probability distributions, and 
the fraction of entangled states at a given time can be 
calculated. However, it is not clear how much entangle- 
ment remains at a given time. This is especially impor- 
tant in the AD case, because the entanglement of most 
states (64.27%) decay asymptotically. In the next three 
sections, we consider the dynamics of the concurrence 
probability distribution for each case. 



We start by investigating the case in which each qubit 
is coupled to a depolarizing channel. As time goes by, 
the depolarizing channels induce decoherence on the two 
qubits and their initial entanglement is lost. The concur- 
rence at time g. 



C^^\q;Co) = max jo, Co (1 - qf - \q{2 - g)| 



(6) 



depends only on (time and) the concurrence of the ini- 
tial state Cq. The loss of entanglement in the system is 
uniform, and initial maximally entangled states remain 
as the states with the highest concurrence for all times, 
Cm [q] = 1 — f 9(2 — (?) • For an ensemble of uniformly dis- 
tributed initial pure states, the concurrence probability 
density at a time q 



P^''HC;q)^3 



C+y{2-q) 



1 - 



g + ^g(2-g) 

(i-qr 



follows from Eq. dH) and p{C; 0) = 3CVl-C^. In the 
inset of Fig. [2] we draw p^{C;q) for different values of 
time, q. 

In Fig. [21 entanglement degradation is apparent. In- 
deed, due to the loss of entanglement, the whole dis- 
tribution moves towards zero and the average concur- 
rence Cq = J^" dC C p^^\C;q) decreases. At the same 
time entanglement concentrates around its mean value 
because the standard deviation (C2 — Cq )^/2 falls as q 
increases. Finally, entanglement is lost: the probability 
to obtain an entangled state j^" dC p^^\C]q) declines 
with time. 

Entanglement degradation, concentration and loss are 
displayed in all the examples considered in this work. 
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Dimensionless Reparameterized Time q 

FIG. 2. (Color online) Identical but independent depolarizing 
channels. The solid blue lines represent the average concur- 
rence plus a standard deviation (darkest line), the average 
concurrence, and the average concurrence minus a standard 
deviation (lightest line); maximum concurrence is plotted in 
red (dashed line), and the value of concurrence at which the 
distribution peaks is plotted in green (dotted line). Inset: 
From right to left, concurrence distributions for q = 0,0.2 
and 0.4. 

and seem typical for open markovian systems, whose 
asymptotic states belong to S and dS. This asymp- 
totic behavior originates also entanglement concentration 
due to the shrinking average distance between the states 
that remain entangled. Entanglement is lost when states 
cross the border between entangled and separable states. 
Thtis, the concurrence probability distribution provides 
a complete picture of the process of entanglement de- 
cay. However, the most salient features of this process 
are contained in the evolution of the concurrence maxi- 
mum, its average, its deviation and the separation-time 
probability distribution. 

V. AMPLITUDE-DAMPING CHANNELS. 

In our second example each qubit experiences dissi- 
pation under the action of an amplitude-damping chan- 
nel. For an initial two-qubit pure state, the concurrence 
evolves in time as [l^l 

^^^) ^ j^^ax {0, (1 - q) (Co - 2|VJn|'g) } . 

(7) 

We notice that the set of initially entangled states sep- 
arates into two classes: those states for which Co < 
2|?/'iip, becoming separable at a finite time; while all 
the others states, for which Co > 2|-0iip, reaching sepa- 
rability only asymptotically. Hence, contrary to the pre- 
vious case, the loss of concurrence is not uniform: when 
comparing states with the same initial concurrence, those 
with smaller are more robust. At a given time q, the 
maximum concurrence Cm = \ — q corresponds to initial 
states of the form (|01) + exp(i0) |10))/-\/2, with arbi- 



trary phase (p. Since C^'^^' {q) depends on both the ini- 
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Dimensionless Reparameterized Time q 



FIG. 3. (Color online) Identical but independent amplitude- 
damping channels. The solid blue lines represent the aver- 
age concurrence plus a standard deviation (darkest line), the 
average concurrence, and the average concurrence minus a 
standard deviation (lightest line); maximum concurrence is 
plotted in red (dashed line), and the value of concurrence 
at which the distribution peaks is plotted in green (dotted 
line). Inset: From right to left, concurrence distributions for 
q = 0,0.4,0.8 and 0.9. 

tial concurrence and the state component in the |11) di- 
rection, the joint probability distribution Co) 
over all two-qubit pure states is calculated first in order 
to evaluate the concurrence distribution p(^°)(C). This 
joint probability either vanish or is equal to 3 Co log((l 4- 
Vl - C^)/z), where z = max(I - - C^f, 2 jiAiip) for 
2|^iiP < 1 -I- ■\/ 1 — Cq . The concurrence probability dis- 
tributionp(AD)(^.g) ^ /d|7/.npdCo<5(/)Pj(|ViiP,Co) 
(where f = C — C^-'^^^q; Co, V'li)) can be calculated an- 
alytically, but the expression is cumbersome and not il- 
luminating. 

The concurrence probability distribution is depicted in 
the inset of Figure [3l for some values of time q. As in the 
previous case the degradation of concurrence (the disen- 
tanglement) is more important for smaller (larger) times. 
The separation between these processes is more conspicu- 
ous here because most states disentangle asymptotically. 
The asymptotic massive disentanglement is not evident 
in the graphic of average concurrence. 



VI. PHASE-DAMPING CHANNELS. 

In our last example each qubit evolves under the action 
of a dephasing channel, which washes away the coher- 
ences in the representation. The concurrence at time 
q of an initial pure two-qubit state under the influence of 
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identical independent phase-damping channels 

{q; Co , r, s) = max 1 0, ,s + ^ + {1 - q) | 

(8) 

with q = 2q-q^, s + r = 2|?Aoo'0ii| and s-r = 2|'0oi'0io|, 
allows us to show that ESD is generic: only a zero- 
measure set of entangled initial states do not disen- 
tangle in finite time, those with at least one vanishing 
component t/i^. The initial states (|01) + e"'' \10)) / V2 
and (|00) -I- e"^ |f f))/-\/2 have the maximum concur- 
rence at time q, (f — q)^. The probability distribution 
of concurrence p^^^^C;q) = J dO ds dr Pj{r, s)5{Cq — 
C^^^\q;Co,r, s))/2tt, depends on the joint probability 
Pj (s,r), defined before. After recasting p(^°)(C; g) as a 
2d integral with complicated integrands and integration 
regions, it is evaluated numerically. The probability dis- 
tribution of concurrence for several values of time q is 
shown in the inset of Fig. 2] 
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FIG. 4. (Color online) Identical but independent phase- 
damping channels. The solid blue lines represent the aver- 
age concurrence plus a standard deviation (darkest line), the 
average concurrence, and the average concurrence minus a 
standard deviation (lightest line); maximum concurrence is 
plotted in red (dashed line). Inset: From right to left, con- 
currence distributions for q = 0, 0.4, 0.57 and 0.8. 



currence show that the disentanglement process occurs 
first for case D (identical but independent depolarization 
channels), then for case PD (phase-damping channels), 
and finally for case AD (amplitude-damping channels), 
where about 2% of the states decay asymptotically. This 
separation order is expected, because it depends on the 
behavior of the set of almost pure states. 
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FIG. 5. (Color online) Disentanglement-time probability dis- 
tributions and average concurrence (Inset): solid line (black) 
for local depolarizing channels (D); dashed line (blue) for lo- 
cal amplitude damping channels (AD); and dotted line (red) 
for local phase damping channels (PD). The delta function at 
q = for the ESD time distribution is common to all cases. 

The behavior of cases D and PD, is quite similar for 
mixed states, up to their different time scales. In fact, 
not only the concurrence averages and the distributions 
of the disentanglement times are similar {p^^^\qs) ~ 

p^^'>{qs/a)/a and C^^°^ w C^°^, where a « 1.56), 
but also the corresponding concurrence distributions are 
alike. As shown in Fig. [6] the probability distributions of 
concurrence satisfy p(^°)(C; g) « p(°^(C; g/a). For early 
times, an appropiate time scaling shows that the concur- 
rence probability distribution of case AD is also similar 
to the other two. For longer times, the relative weight of 
small values of concurrence is higher in this case. 



VII. MIXED STATES. 

Diferent experimental conditions call for different sets 
of states to be considered, in particular mixed states. 
Albeit no natural uniform distribution exists for them, 
Hilbert-Schmidt and Bures measures are commonly used. 
Analytical results can not be obtained along the lines 
presented in the main text, but it is possible to per- 
form numerical simulations. Employing the first mea- 
sure, around 24% of the states are separable (before in- 
teraction with the channels), as signaled by a delta func- 
tion contribution at time g = in Fig. [5l The probability 
distribution of the separation times and the average con- 



VIII. A SINGLE NOISY CHANNEL. 

In order to have a wider spectrum of possible behaviors 
we also consider the case in which only one of the qubits 
is affected by a dissipative channel (D, AD, PD), and the 
initial states are pure. The concurrence at time q turns 
out to be C(g) = max(0, x(g)Co) where x^^\q) = 1 — 
3q/2, x'-^^^q) = yr^andx^^^n?) = l-?- If we have 
a single depolarizing channel, the process of concurrence 
degradation occurs up to the reparameterized time q = 
2/3, and then all states undergo entanglement sudden 
death, but in the other two cases concurrence decreases 
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Concurrence C 

FIG. 6. (Color online) Concurrence distributions: solid line 
(black) for local depolarizing channels (D) and times q — 
0,0.15,0.31; dashed line (blue) for local amplitude damping 
channels (AD) and times g = 0, 0.2, 0.55; and dotted line 
(red) for local phase damping channels (PD) and times q — 
0,0.1,0.2. The rightmost distribution (solid) is common to 
all cases. 

steadily to zero in an asymptotic way, and at a slower 
pace in the case of the amplitude-damping channel. 

IX. NON-AUTONOMOUS SYSTEMS 

We also consider non-autonomous systems [sj], whose 
equation of motion is of the form p = exp(—jt)C{p). If 
the Kraus operators of the usual dissipative dynamics 
p = C{p) are Ei{q), the Kraus operators of the non- 
autonomous version are Ei{q/"f). That is, the dynamics 
at infinite time is given byp(l/7): if7 > l,allofthc cases 
considered in this paper would suffer ESD, but if 7 < 1. in 
the cases of single or double AD and PD channels, there 
would be a non-zero-measure set of asymptotic entangled 
states. Analogously, we would have asymptotic entangled 
states for the single depolarizing channel if 7 < 2/3, and 
for two identical depolarizing channels if 7 < 1 — l/\/3- 

X. TIME-DEPENDENT MARKOVIAN AND 
NON-MARKOVIAN CHANNELS 

As demonstrated in [27| an infinite number of open 
quantum systems can be described by the same Kraus 
representation used here to exemplify our results in a 
Markovian, constant decay-rate setting. All of our results 
can be mapped to any of these time-dependent Marko- 
vian or non-Markovian channels, as we show below. In 
order to fix ideas, we describe the interaction of a single 
qubit interacting with its environment by the Hamilto- 
nian 

H = hujocr^a^ + hwkO^f^ak + {o+A + (t_j4^), (9) 

k 



with A = '^i^.figkak, where Qk are coupling constants 
and h is Planck's constant. We see the qubit do not 
need to be (strictly) a quantum memory, due to the first 
term of the Hamiltonian. The environment is modelled 
as a collection of quantum oscillators of frequency ujk, 
whose annihilation and creation a\. operators satisfy 
the usual boson commutation relations. The transition 
frequency of the qubit is denoted by wq, and a± are the 
qubit raising and lowering operators. The initial state 
of the total system, |-0) (g) Ylk is separable, with the 
oscillators in their ground state. For positive times, the 
system's state is 

\^it)) = a Iff, 0) + /3 (^coit) |e, 0> + E c„(t) \g, 1„)^ , 

where |0) = HJO),, |1„) = (n^^jO)^) ® |1)„, and 

co{t) and c„(t) are time-dependent coefficients such that 
co(0) = 1 and c„(0) = 0. The Schrodinger equation 
of motion provides equations for the time derivatives of 
these cocfhcients. Formally solving the equations for the 
coefficients c„(t), and replacing back into the equation 
for co(t), we obtain the integro-differential equation 

co{t) + zc^oco(t) + fdrj^ |5«l'e-'""(*-")co(T) = 0. 

In the continuum limit, the environment's correla- 
tion function |(7„pe~'""'*~'^'' becomes the Fourier 
transform of the spectral density J(w). that is, 

J(^)e--(*-). Choosing J{u;) = ^^j.^^^, 

where A^^ is the environment's correlation time and 70 is 
the system's relaxation rate in the Born-Markov approx- 
imation, allows for an analytic solution for the coefficient 
co(i), dl 

co{t) = e~^-'"°* ( cos(m) + A sm(17t) ) , (10) 



where ft — ^ A(27o — A)/2. Tracing out the environ- 
mental degrees of freedom we obtain the qubit dynam- 
ics, which corresponds to an amplitude-damping (AD) 
channel with Kraus operators Ei(t) = e~*'^''*l^^^^l_Ei, 
and i?o(t) = |0)(G| + ^\ - q{i)\\){\\ and Ex[t) ^ 
A/g(t) |0)(1|. The time-dependent q(t) = 1- |co(t)P goes 
from zero at the initial time to one at large times, but 
it has an infinite number of intervals where it decreases 
from one down to a local minimum, and back to one. See 
the dashed line in Figure [71 where 70 was chosen to be 
equal to 4A. Since we know the behavior of entangle- 
ment as a function of q, we also know it in function of t. 
In particular we see that, the intervals of decreasing g(t) 
correspond to negative decay rates, increase of average 
entanglement and entanglement dispersion, and sudden 
birth of entanglement. If 270 < A, contrary to what has 
been assumed here, the dynamics is Markovian, but the 
decay rate is not constant. 



7 




because it assesses 



2 4 6 

Physical Time t (in units of 1/A) 

FIG. 7. (Color online) Typical forms of the parameterized di- 
mensionless time q{t). From left to right: dashed line (blue) 
for a non-Markovian channel; solid thin line (red) for time- 
dependent Markovian channel; solid thick line (black) for a 
non-Markovian periodic channel; and dotted line (black) for 
a non-Markovian time-dependent channel with incomplete de- 
coherence for long times. Time is measured in units of 1/A, 
where A is a frequency, as described in the text. 



It is possible to formulate time-dependent Marko- 
vian and non-Markovian models for all the noisy chan- 
nels considered here (see [l^l). All the details of these 
non-constant decay rate channels are completely cap- 
tured in q{t). Some typical dependencies, for a dephas- 
ing model where q{t) = 1 — exp(— r(i)) and r(t) = 
/ du J{uj) coth{fuv/2kBT), are depicted in Fig. 

[71 Here, Bolztmann's constant is denoted by fcs and 
the temperature of the environment by T. The thick 
solid line corresponds to a single oscillator, where q{t) 
oscillates periodically and, hence does not have a limit 
for long times. The thin solid line and the dotted 
line correspond to an ohmic environment, for which 
J(a;) = (w/A) exp(— cj/A), for finite temperature and zero 
temperature, respectively. While both exemplify time- 
dependent Markovian processes, the latter presents in- 
complete decoherence, for which q{t) does not approach 
one for long times. Note that non-Markovian channels 
present a more varied behavior than the non-autonomous 
channels considered in the previous section. 



XI. CONCLUSIONS AND OUTLOOK. 

This research shows that concurrence probability 
and separation-time distributions constitute a new, 
global approach to the decay of quantum entangle- 



ment, exemplified here in the case of two-qubit quan- 
tum memories under identical local decoherence mecha- 
nisms. The latter distribution follows from the former 

t=ts 

the weight of the separable states. We have observed 
that entanglement not only disappears suddenly but also 
degrades and converges to its mean in a characteristic 
way for different decoherence channels, when the initial 
states are assumed to be pure and uniformly distributed. 
Despite this convergence to the mean, the maximum con- 
currence is usually much higher than the average one. 
Thus, whenever the initial state can be chosen, e.g. in 
teleportation protocols as a quantum resource, using the 
most robust states offers a great advantage. 

Several variations of the main problem considered in 
this paper were analyzed. First, if the initial states are 
mixed, the global concurrence evolution is similar in all 
cases, except for a different time-scale, up to times where 
the effect of the most entangled states becomes impor- 
tant. The posibility of a universal behavior for initial 
mixed states must be further explored. Second, if the 
initial states are pure, but only one qubit suffers deco- 
herence, no state separates at a finite time or all of them 
separate at the same time. Finally, if the evolution of the 
qubits is non-Markovian or non-autonomous (3l| , our re- 
sults must be changed. For example, it may happen that, 
for very long times q{t) does not reach 1. Depending on 
the value reached, (significant) entanglement might re- 
main for long times. In non-Markovian processes q(t) 
might also oscillate; hence, as times goes by, entangle- 
ment not only dies suddely (degrades, concentrates) but 
also appears suddely (becomes enhanced, disperses). 

The derivation of entanglement statistics for interact- 
ing multiqubit systems under the action of collective or 
individual noisy channels, and the investigation of the 
consequences in quantum information processing, present 
an interesting challenge for future research. 
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